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We study several cosmological models with Bianchi VIo&III symmetries under the self-similar 
approach. We find new solutions for the "classical" perfect fluid model as well as for the vacuum 
model although they are really restrictive for the equation of state. We also study a perfect fluid 
model with time-varying constants, G and A. As in other models studied we find that the behaviour 
of G and A are related. If G behaves as a growing time function then A is a positive decreasing time 
function but if G is decreasing then Ao is negative. We end by studying a massive cosmic string 
model, putting special emphasis on calculating the numerical values of the equations of state. We 
show that there is no SS solution for a string model with time- varying constants. 
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I. INTRODUCTION. 

The current observations of the large scale cosmic mi- 
crowave background suggest us that our physical universe 
is expanding, isotropic and homogeneous models with 
a positive cosmological constant. The analysis of CMB 
fluctuations could may confirm this picture. But other 
analysis reveal some inconsistency. Analysis of WMAP 
data sets show us that the universe could have a preferred 
direction. For this reason Bianchi models are important 
in the study of anisotropies. 

The study of SS models is quite important since a 
large class of orthogonal spatially homogeneous models 
are asymptotically self-similar at the initial singularity 
and are approximated by exact perfect fluid or vacuum 
self-similar power law models. Exact self-similar power- 
law models can also approximate general Bianchi models 
at intermediate stages of their evolution. This last point 
is of particular importance in relating Bianchi models to 
the real Universe. At the same time, self-similar solu- 
tions can describe the behaviour of Bianchi models at 
late times i.e. as t — > oo (see [l|).Of particular interest 
is determining the exact value of the equations of state 
(see @). The geometry and physics at different points 
on an integral curve of a homothctic vector field (HVF) 
differ only by a change in the overall length scale and 
in particular any dimcnsionless scalar will be constant 
along the integral curves. In this sense the existence of a 
HVF is a weaker condition than the existence of an KVF 
since the geometry and physics are completely unchanged 
along the integral curves of a Killing vector field (KVF) . 
However, the existence of a non-trivial HVF leads to re- 
strictions on the equations of state. We shall put special 
emphasis on this fact. 

In modern cosmological theories, the cosmological con- 
stant remains a focal point of interest (see for 
reviews of the problem). A wide range of observations 
now compellingly suggest that the universe possesses a 
non-zero cosmological constant. Some of the recent dis- 
cussions on the cosmological constant "problem" and 
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on cosmology with a time-varying cosmological constant 
point out that in the absence of any interaction with 
matter or radiation, the cosmological constant remains 
a "constant" . However, in the presence of interactions 
with matter or radiation, a solution of Einstein equa- 
tions and the assumed equation of covariant conserva- 
tion of stress-energy with a time- varying A can be found. 
This entails that energy has to be conserved by a de- 
crease in the energy density of the vacuum component 
followed by a corresponding increase in the energy den- 
sity of matter or radiation Recent observations strongly 
favour a significant and a positive value of A with mag- 
nitude A(Gh/c 3 ) w 10~ 123 . These observations suggest 
an accelerating expansion of the universe, q < 0. 

Therefore the paper is organized as follows. In section 
II we outline the metrics as well as the main geometrical 
ingredients of the models, as for example, all the cur- 
vature invariants. In section III we start by reviewing 
some basic notions on self-similarity (SS). We calculate 
the HVF for the metric and obtain the restrictions for 
the scale factors. Section IV is devoted to revising the 
"classical" solution for a perfect fluid (A = 0). Neverthe- 
less we find a new solution in each of the studied cases 
i.e. Bianchi types III and VI. In section V we shall study 
the vacuum solution, as in the above section we obtain 
solutions for Bianchi types III and VI. Once we have fin- 
ished with the classical solutions we go next to the study 
of SS solutions in different contexts. For this purpose, 
in section VI we study a perfect fluid model with time- 
varying constants, G and A. As we have mentioned above 
we shall pay special attention to calculating the numeri- 
cal value of the equation of state. In section VII we study 
a massive cosmic string model while in section VIII we 
study its generalization i.e. by considering a model which 
allows time-varying constants. In the last section we end 
by summarizing the main results. 



II. THE METRICS AND THE GEOMETRIC 
INGREDIENTS. 

Throughout the paper M will denote the usual 
smooth (connected, Hausdorff, 4-dimensional) spacetime 
manifold with smooth Lorentz metric g of signature 



2 



(—,+,+,+) (see for example [Ti|). Thus M is paracom- 
pact. A comma, semi-colon and the symbol C denote the 
usual partial, covariant and Lie derivative, respectively, 
the covariant derivative being with respect to the Levi- 
Civita connection on M derived from g. The associated 
Ricci and stress-energy tensors will be denoted in com- 
ponent form by R%j{= R c j C d) and T^- respectively. We 
shall use a system of units where c = 1. 

A Bianchi VI space-time is a spatially homogeneous 
space-time which admits a group of isometries G 3 , acting 
on spacelike hypersurfaces, generated by the spacelike 
KVs, 

£1 = d x , £ 2 = d y , £3 = mxd x - nydy + d z , (1) 

and therefore C^ 3 = m, C| 3 = — n. In synchronous co- 
ordinates the metric is: 



df + a\t)e- lmz dx A + b l {t)e Znz dy z + d\t)dz l 



ds z 

(2) 

where the metric functions a(t),b(t),d(t), are functions 
of the time co-ordinate only. 

As it is observed, the metric © collapses to the fol- 
lowing cases: 

1. if m = —n, then the metric collapses to Bianchi 
type V, 

ds 2 = -dt 2 + a 2 (t)e- 2mz dx 2 + b 2 (t)e- 2mz dy 2 + d 2 (t)dz 2 , 

(3) 

with 

£1 = d x , £ 2 = d y , £3 = mxd x + mydy + d z , (4) 
and therefore C\ z — m = Cf 3 . 

2. if n = 0, then the metric collapses to Bianchi type 
III, 

ds 2 = -dt 2 + a 2 (t)e- 2mz dx 2 + b 2 (t)dy 2 + d 2 (t)dz 2 , (5) 
with the following KVF (Killings vector fields) 

£1 = d x , £2 = d y , £ 3 = mxd x + d z , (6) 
and therefore Cf 3 = m. 

3. If m = n = 0, then the model collapses to Bianchi 
I model, 

ds 2 = -dt 2 + a 2 (t)dx 2 + b 2 {t)dy 2 + d 2 (t)dz 2 , (7) 

with the following KVF (Killings vector fields). 

ti = d x , & = d y , & = d z , (8) 

and therefore C£ = 0. 

Once we have defined the metric and we know which 
are its killing vectors, then we calculate the four velocity. 



It must verify the equation, C^Ui = 0, so we may define 
the four velocity as follows: 

u l = (1,0,0,0), (9) 

in such a way that it is verified, g(u\ u l ) = — 1. 
From the definition of the 4— velocity we find that: 
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H 



q= dt\H ] - 1 ' 



(10) 



We shall take into account the Einstein's field equa- 
tions (FE) written in the following form: 



Ri 



1 



Rgij = 8irGTij — Ag 



(11) 



where, Ty, is the energy-momentum tensor. 

We also study the curvature behaviour of the different 
solutions (see for example [H, 0, 03 and (3)- The 
studied curvature quantities are the following ones. Rice 
Scalar, I — R\, Krestchmann scalar, l\ := iJyjyiJ t3 ' fc ', 
the full contraction of the Ricci tensor, I2 := RijR 1 ^. 
The non-zero components of the Weyl tensor. The Weyl 
scalar, 7 3 = C abcd C abcd = h - 2I 2 + as well as the 
electric scalar I4 — EijE^ , (see [3) and the magnetic 
scalar ^5 = if,j , of the Weyl tensor. The Weyl param- 
eter (see W.C. Lim et al [lg|) which is a dimensionless 
measure of the Weyl curvature tensor 



W 2 = 



H 4 



6H 4 V y 



//-.-//'") 



24H 4 ' 



(12) 



W can be regarded as describing the intrinsic anisotropy 
in the gravitational field. And to end, we shall calcu- 
late the gravitational entropy. From a thermodynamic 
point of view there is every indication that the entropy 
of the universe is increasing. Increasing gravitational 
entropy would naturally be reflected by increasing local 
anisotropy, and the Weyl tensor reflects this. One sug- 
gestion in this connection was Penrose's formulation of 
what is called the Weyl curvature conjecture (WCC) [l7l |. 
The hypothesis is motivated by the need for a low en- 
tropy constraint on the initial state of the universe when 
the matter content was in thermal equilibrium. Penrose 
has argued that the low entropy constraint follows from 
the existence of the second law of thermodynamics, and 
that the low entropy in the gravitational field is tied to 
constraints on the Weyl curvature. Wainwright and An- 
derson [lll | express this conjecture in terms of the ratio 
of the Weyl and the Ricci curvature invariants, 



P 2 



(13) 
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The physical content of the conjecture is that the ini- 
tial state of the universe is homogeneous and isotropic. 
As pointed out by Rothman and Anninos (l9l . I20I ] the en- 
tities P 2 and I3 are "local" entities in contrast to what 
we usually think of entropy. Gr0n and Hervik ( [ijl, 0] ) 
have introduced a non-local entity which shows a more 
promising behaviour concerning the WCC. This entity is 
also constructed in terms of the Weyl tensor, and it has 
therefore a direct connection with the Weyl curvature 
tensor but in a "non-local form". 

For SS spacetimes, Pelavas and Lake ([12]) have 
pointed out the idea that eq. (flU)) is not an acceptable 
candidate for gravitational entropy along the homothctic 
trajectories of any self-similar spacetime. Nor indeed is 
any "dimensionless" scalar. It is showed that the Lie 
derivative of any " dimensionless" scalar along a homoth- 
etic vector field (HVF) is zero, and concluded that such 
functions are not acceptable candidates for the gravita- 
tional entropy. Nevertheless [23j , since self-similar space- 
times represent asymptotic equilibrium states (since they 
describe the asymptotic properties of more general mod- 
els), and the result P 2 — const. , is perhaps consistent 
with this interpretation since the entropy does not change 
in these equilibrium models, and perhaps consequently 
supports the idea that P 2 represents a " gravitational en- 
tropy". As we shall show W 2 and P 2 will be constant 
along homothetic trajectories, since all the dimension- 
less quantities remain constant along timclikc homothetic 
trajectories. 

III. THE SELF-SIMILAR SOLUTION. 

In general relativity, the term self-similarity can be 
used in two ways. One is for the properties of space- 
times, the other is for the properties of matter fields. 
These are not equivalent in general. The self-similarity in 
general relativity was defined for the first time by Cahill 
and Taub [24| , Eardley [25[ (see for general reviews [2t| - 
(28|). Self-similarity is defined by the existence of a ho- 
mothetic vector V in the spacetime, which satisfies 

C v g tj = 2ag tj , (14) 

where g^ is the metric tensor, Cy denotes Lie differenti- 
ation along V and a is a constant. This is a special type 
of conformal Killing vectors. This self-similarity is called 
homothety. If a ^ 0, then it can be set to be unity by a 
constant rescaling of V. If a = 0, i.e. CyQij = 0, then V 
is a Killing vector. 

Homothety is a purely geometric property of space- 
time so that the physical quantity does not necessarily 
exhibit self-similarity such as LyZ = kZ, where k is a 
constant and Z is, for example, the pressure, the energy 
density and so on. From equation (fT4"|) it follows that 
CyR 1 jki — 0, and hence CyRij — 0,and LyGij = 0. A 
vector field V that satisfies the above equations is called 
a curvature collineation, a Ricci collineation and a mat- 
ter collineation, respectively. It is noted that such equa- 



tions do not necessarily mean that V is a homothetic vec- 
tor. We consider the Einstein equations Gij — 8wGTij, 
where Ty is the energy-momentum tensor. If the space- 
time is homothetic, the energy-momentum tensor of the 
matter fields must satisfy CyTij — 0. For a perfect 
fluid case, the energy-momentum tensor takes the form 
Tij — (p+p)uiUj+pgi 3 where p and p are the pressure and 
the energy density, respectively. Then, equations UJJ re- 
sult in 

Cyu 1 — —au l , Cyp = —Zap, CvP = —2ap. (15) 

As shown above, for a perfect fluid, the self-similarity 
of the spacetime and that of the physical quantity co- 
incide. However, this fact does not necessarily hold for 
more general matter fields. Thus the self-similar variables 
can be determined from dimensional considerations in the 
case of homothety. Therefore, we can conclude homoth- 
ety as the general relativistic analogue of complete simi- 
larity. From the constraints (|15p , we can show that if we 
consider the barotropic equation of state, i.e., p = f(p), 
then the equation of state must have the form p — up, 
where w is a constant. In this paper, we would like to try 
to show how taking into account this class of hypothesis 
one is able to find exact solutions to the field equations 
within the framework of the time varying constants. 

From equation (fT4"|) ,we find the following homothctic 
vector field for the BVLj metric @: 

V = td t + (\ - t~^j xd x + (\ - t^j yd y + (\ - t~ J zd z , 

(16) 

with the following constrains for the scale factors: 

a(t) = a t ai , b(t)=b t a2 , d{t) = d t, (17) 

where 01,02 6 t and we set do = 1, therefore the result- 
ing homothetic vector field is: 

V = td t + (l-a 1 )xd x + {l-a 2 )yd y , (18) 

so the metric © collapses to the following one: 

ds 2 = -dt 2 + a 2 {t)e- 2mz dx 2 + b 2 (t)e 2nz dy 2 +t 2 dz 2 . (19) 



IV. PERFECT FLUID SOLUTIONS. 

The energy-momentum tensor for a perfect fluid 
model, T^, is defined as follows: 

= (p + p)uiUj + pg i3 , (20) 

with the usual equation of state p = top, u> £ R. The 
resulting FE for the metric with a perfect fluid model 



4 



yield (with A = 0): 

a'b' a! d! d'b' K 



a b ad d b d? 



8irGp, 



, a' d'\ (d! b'\ 
m[ ^--d) +n \-d~-b ]= ^ 



(21) 
(22) 



b" d" d'b' n 2 , s 



d" 



a' d' m 



d a ad d? 
b" a" a'b' ran 
b a a b d 2 



= -SuGujp, (24) 



really restrictive since it only has physical meaning if 
w < 0. 

For this solution we have obtained the follow- 
ing behaviour for the main curvature quantities. 
We start by calculating the curvature invariants, 
(Iif i=0 , they yield: /„ = 2 (a| + a 2 - m 2 + l) t~ 2 , 
h = 4 (af - 2a| + 3a 2 , + to 4 - 2m 2 + l) i~ 4 , and I 2 = 
2 [a^ + 2a\ + 2a 2 + m A -2m 2 + 1 - 2m 2 a 2 ) t~ i . As it 
is observed we have obtained a singular solution. 
The non-zero components of the Weyl tensor are: 
C txtx = -(l/6)Ke- 2m * : C tyty = (l/3)0 2 ^- x ), 



SirGtup, (25) C tztz = -(l/6)/C, = (l/6)/C e - 2m ^ 2a2 , 



, , a' 6' d' , 

P ' + P (l + .)[- + - + -\=0, 



(26) 



2 ) 



Therefore, taking into account the restrictions for the 
scale factors given by eq. (jTTJ) , the solution for a perfect 
fluid model, eqs. (|2HI26p , is the following one. 

Form (|2"6"|) we get: 

p = Po t-( 1 +^)( a i+^+ 1 ) = Po t-y, 

with (1 + u>) (ai + a 2 + 1) = (1 + uj) a = 7. Now, taking 
into account eq. pTj) we get aia 2 +a 1 + a 2 + K = SirGpo, 
and (1 + uj) (ai + a 2 + 1) = 2, and hence (|2HI26p collapse 
to the following algebraic system of equations: 

m(ai-l)+n(l-a 2 ) = 0, (27) 

a 2 (a 2 - 1) + a 2 - n 2 = -Aw, (28) 

01 (ai - 1) + ai - to 2 = -Ah), (29) 

^2 («2 — 1) + ai (ai — 1) + aia 2 + ran = — Aj, (30) 

(1 + w) (01 + oa + 1) = 2, (31) 

with 

.4 



Po 



A = aia 2 + ai + a 2 + K, 



8ttG' 

obtaining the following set of solutions. 

A. Bianchi type III. 

Our first solution is the following one: 
2w 



ai = 1, a 2 



uj + 1 



V-3uj 2 + 2uj + 1 

m = , n = 0, 

uj + 1 



(32) 

note that the solution only has sense if uj E (— 5,1), 
and therefore we obtain the following values: a% = 1, 
a 2 G (0, 1) , iff w < 0, m G (0, 1) ,and (1 + w) (2 + a 2 ) = 
7 = 2, so p = pot~ 2 , and therefore the metric collapses 
to the following one 



-dt 2 + t 2 e~ 2mz dx 2 + t 2a2 dy 2 + t 2 dz 2 , (33) 



ds 2 



which stand for a Bianchi III metric. Note that 
admits the KVF given by eq. ©• This solution looks 



-(l/3)/Ce- 2mz i 2 , and C yzyz = (l/6)0 2a2 , where 
K. = (to — 1 + a 2 ) (m + 1 — a 2 ) . Therefore the Weyl in- 
variant yields, I3 = (4/3)/C 2 i~ 4 . The electric invariant 
yields: I4 — (l/6)/C 2 i~ 4 , while the magnetic invari- 
ant, 1$, vanishes. The Weyl parameter yields, W 2 = 
/C 2 /36(2 + a 2 ) 4 = const., while the gravitational entropy 
yields, P 2 — const., as it is expected. 



B. Bianchi type VIo, solution I. 

The classical solution given by Collins (29|, Hsu et al 
(see [32|), i.e. 



a\ = a 2 



1-uj 
2(uj + 1) 



V-3uj 2 + 2uj + 1 

m = n= — — , (34) 

2 (uj + 1) 



with we (— |, l) , and therefore the metric collapses to 
the following one 



ds 2 = -dt 2 + i 



2ai 



-2rnz 



dx 2 +t 2ai e 2mz dy 2 +t 2 dz 2 , (35) 



finding in this way that, a 2 G (0, 1) , m G (0, and 
(1 + uj) (1 + 2a 2 ) = 7 = 2, so p — pot~ 2 - In the way we 
find that 



„ 1 + 2a 2 

H = , q 



-2a 2 



l + 2a 2 



<0, 



(a 2 - l) 2 
3t 2 



For this solution we have obtained the follow- 
ing behaviour for the main curvature quantities. 
We start by calculating the curvature invariants, 

(Ji)L>> the y y ield: J o = 2 (3a 2 - to 2 ) t- 2 , h = 

4 (3a 4 - Aa\ + Aa\ - 2a 2 TO 2 + 3to 4 + 4toV - 4to 2 ) i~ 4 , 



and I 2 = 2 (3a 4 - 2a{ + 2a{ 



2m 2 a\) t 4 . As it is 



observed we have obtained a singular solution, as it is 
expect for a SS solution. The non-zero components of the 
Weyl tensor are: C txtx = {l/3)m 2 t 2 ^-^ e~ 2mz , C txxz 



m(l - a^t 2 ^- 1 



-2mz 



tyty 



(l/3)TO 2 t 2 ( ai ~ 1 )e 2mz , 



TO(a 1 -l)i 2a i- 1 e 2mz , C tztz = -(2/3)m-, 



C xyxy = (2/3)m 2 t 2 ( 2 ^- 1 ), C x 



■(1/3)1; 



2a i 



I 2 , 
-2m z 



and C yzyz = (l/3)t 2ai e 2mz . Therefore the Weyl in- 
variant yields, I 3 = (16/3) (-3a 2 + 6ai - 3 + to 2 ) i" 4 . 
The electric invariant yields: I4 — (2/3)TO 4 t -4 , 
while the magnetic invariant, ^5, behaves as: 
Is = 2m 2 (ai — I) 2 t~ 4 . The Weyl parameter yields, 



■5 



i 2 (3a 2 



6ai + 3 + m 2 )) /9(1 + 2ai) 4 = const., 



while the gravitational entropy yields, P 2 
is expected. 



const., as it 



C. Bianchi type VIo, solution II. 

We find a new solution, where 

1 — (l 2 (U) + 1) — CjJ 

ai = — , a 2 = o 2 , 



uj + 1 
6 (q 2 (a; + 1) + 2u) 

(wTT) ' 



m = 6(l-o 2 ), (36) 



1 - w 
3wTT' 



and at a first look we may say that to e (— 1,1], and 
therefore be [0, oo) . 

The metric collapses to the following one 

ds 2 = -dt 2 + t 2ai e~ 2mz dx 2 + t 2a2 e 2nz dy 2 + t 2 dz 2 , (37) 

finding in this way that, ai = a\ (a 2 ,u>) , so depending 
on the different values of u> we get: 



UJ 


ai 


a 2 


n 


ra 


1 


-a 2 


a 2 








1/3 


\ -02 


(0,1/2) 


^ (2«2 + 1) 


^(l-a 2 ) 





1 — a 2 


(0,1) 


»2 


1 — a 2 


- ("l/3) + 


2-a 2 


(0,2) 


n 


m 



Note that the case u) = 1 looks nonphysical since 
the scale factor a x < 0, for this reason we shall not 
take into account this particular case. Therefore the 
solution is only valid if u) G (— |, l) and hence we get 
(1 + lo) (1 + a\ + a 2 ) = 7 = 2, so p = pot~ 2 . In the way 
we find that 



ff = 



1 + oi + a 2 

t : 



Q = 



- (di + Qg) 
1 + ai + a 2 



<0, 



and 



3t 2 



1 — aia 2 — ai — a 2 ) 



For simplicity we calculate the curvature behaviour 
for the case u> = 0. We start by calculating the curva- 
ture invariants, (ii) i=0 , they yield: 7 = 4a 2 (a 2 — 1) t~ 2 , 
/i = 16a 2 (a 2 — l) 2 i~ 4 = 7 2 . As it is observed we have 
obtained a singular solution, as it is expect for a SS so- 
lution. The non-zero components of the Weyl tensor are: 

, J. — 2(12 r' ~, +1— 2a2 /"< +2(02 — 1) /"< 

Ct^tz = const., Cxyxy = const., C^z ~ i 2(1 ~ a2 \ 



2o 2 -l 



f 

and C yzyz <~ t 2 " 2 . Therefore the Weyl invariant yields, 
73 = 32a 2 , (a 2 — 1) /3t 4 . The electric invariant yields: 
h = 2a\ (o 2 - 1) /3t 4 , while the magnetic invariant, Z5, 
I5 = 2a 2 . (a 2 — l) 2 t~ 4 . The Weyl parameter 
= a\ (a 2 — l) 2 /36 = const., while the gravi- 



behaves as: 
yields, W 2 



SS VACUUM SOLUTION. 



In this case, the system to solve is the following one: 



ai<2 2 + ai + a 2 + (ran — m 2 — n 2 ) = 0, 
m (ai — 1) + n (1 — a 2 ) = 0, 
a 2 (a 2 - 1) + a 2 - ra 2 = 0, 
ai (ai — 1) + ai — to 2 = 0, 
a 2 (a 2 — 1) + ai (ai - 1) + a x a 2 + ran = 0, 



(38) 
(39) 
(40) 
(41) 
(42) 



obtaining the following solutions (as well as the trivial 
one). 



1. We find our first solution as 
a 2 = a 2 , n = —a 2 , 



m, 



(43) 
(44) 



a, = - I I t v -2aj+4a 2 + l 

and where the solution only has math sense if 

a 2 E ^1 — ^V6, 1 + ^V6 

the negative solution may have sense in the study of 
singularities. We only consider the positive range 
i.e. a 2 G (0, 2.2247) , and therefore we get, a\ G 
(0,1,35). 

Equivalcntly we find 



a 2 = a 2 , 



a 2 , 



-2a\ + 4a 2 + 1 



-m, 



so the metric collapses to the following one 

ds 2 = -dt 2 + t 2ai e 2aiz dx 2 + t 2a2 e 2a2Z dy 2 + t 2 dz 2 , (45) 

it describes a metric belonging to Bianchi type VI 
(i.e. it admits the KVF given by 



£ 3 = -a x xd x - a 2 d y + d z , 

1 - - -a 2 , and to the best 



tational entropy yields, P = 2/3., as it is expected. 



£l — 9 X , £ 2 — ISy, 

and hence Cj 3 = — a\, C 23 
of our knowledge is new. 

For this solution we have found a really curious 
curvature behaviour. The curvature invariants 
(Ii)i— , vanish. The non-zero components of the 
Weyl tensor are: C txtx = - (1/2) /Ce 2 " 1 ^ 2 ^ 1 " 1 ), 
C txxz = (l/2)JCe 2 ^ z t 2 ^-\ C tyty 
(l/2)!Ce 2a2Z t 2< - a2 - 1 \ Ctyyz - - (1/2) JCe 2a2Z t 2a2 -\ 
C xzxz = -{l/2)JCe 2a ^t 2a \ and C yzyz = 
(l/2)JCe 2a2Z t 2a2 7 where JC = (a a - 1 + a 2 )(oi - a 2 ). 
Therefore the Weyl invariant yields, I 3 = 0. 
The electric invariant yields: I 4 = (l/2)/C 2 t -4 , 
while the magnetic invariant, 1$, behaves as: 
h = (l/2)JC 2 t- 4 . The Weyl parameter yields, 

W 2 = K 2 j ^6 (01 + a 2 + l) 4 ^) = const., while the 

gravitational entropy yields, P 2 = 00, as it is 
expected, since I 2 = 0. 
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2. And the set of solutions: 



«2 


= n 


= 0, 


«i 


= m 


= 1, 




«2 


= n 


= 1, 


a i 


= m 


-o, 


and 


«2 


= n 


-o, 


a i 


= 1, 


m 


= -1, 



so the metric collapses to the following one (see Hsu 
et al. [HJ section 2.6.3) 

ds 2 = -dt 2 + t 2 e 2z dx 2 + dy 2 + t 2 dz 2 . (46) 

This solution describes a Bianchi III metric, admit- 
ting three KVF given by eq. © with rn = — 1. For 
this solution we have found a really pathological 
behaviour. All the curvature quantities vanish. So 
this solution lacks of any physical interest. 

VI. PERFECT FLUID MODEL WITH 
VARIABLE CONSTANTS. 



With all these result we find that the system to solve 
is the following one: 

m(a 1 -l)+n(l-a 2 )=Q, (58) 

a 2 {a 2 - 1) +a 2 - n 2 = A, (59) 

ai (ai - 1) + a x - m 2 = A (60) 

02(^2 — 1) + + Oiiip-i — 1) + win — A, (61) 

where A = A whose solutions are the following 

ones. 



A. Bianchi type VIo. 

As it is observed we obtain a solution that is valid for 
all equation of state uj 1 with 

a\ + 1 

a 2 = a 2 , ai = — , 

a 2 + 1 



The resulting FE for the metric ((2|) with a perfect fluid 
matter model (|20j) and with G and A time-varying are: 



a' 


6' 


a' 




a 




h - 

a 














b" 


5) 


+ 




T 


+ 





d' 



A" 

rf2 



— = 8nGp + A, 



m 



b' 



" l d-b 



= 0, 



~~d 

b 



a 
a 

a 



n 

d 2 
,2 



d'b' 
f ~d~b 
a' d' m" 
ad d? 
a' b' mn 
~a~~b + ~d T 
p(l+cu)H 

A' 



-8irGup + A 

~8nGuop + A 

SttGlop 
0, 

-8ttG' P . 



A. 



(47) 
(48) 

(49) 
(50) 

(51) 

(52) 
(53) 



Now, we shall take into account the obtained SS re- 
strictions for the scale factors given by eq. (fT7|) . 
From eq. (fB"2"]) we get 



P = Pot 7 , 

where 7 = (u + 1) a and a = (a± + a 2 + 1). 
From eq. (|4T|) we obtain: 



A 



Ar 



%TTG P0 t- {u+1)a 



(54) 



(55) 



where A — aia 2 +ai+a 2 + K, and K — {run — m 2 — n 2 ) . 

Now, taking into account eq. (|5"3"|) and eq. ([53)) . algebra 
brings us to obtain 



G = Gnf 1 - 



Go = 



A 



4-Kpo (ii) + 1) OL 

While the cosmological "constant" behaves as: 



= A r 



(56) 



(57) 



1 - 2a| 



2a 2 + 1 



ma 2 
al+T 



So this solution only has physical meaning if a 2 € 
(0, 1) , and therefore we have : a\ £ (0.82, 1) , m G (0, 1) , 
n € (0, 1) , so the metric takes the following form: 



ds 2 



-dt 2 + t 2ai e~ 2mz dx 2 + t 2a2 e 2nz dy 2 + t 2 dz 2 . (62) 



Therefore we have the following behaviour for the 
main quantities:. 7 = (u> + 1) a G (0.6,6), a = 
(ai + a 2 + 1) G (1.82, 3) , and A = a 1 a 2 +a 1 +a 2 +K > 0, 
therefore 

P = Po^ 7 , 7 G (0.6, 6), 

is a positive decreasing time function. 
G behaves as: 

G = G f<- 2 , 

so it will be a growing function if 7 > 2, constant if 7 = 2 
and a decreasing time function if 7 < 2. 
The cosmological constant behaves as 



A = ^l--j^, 

hence its sign will depend on the value of 7. If 7 < 2 then 
we get a negative A, it will behave as a true constant if 
7 = 2, and we find that it is a positive time decreasing 
function if 7 > 2. Note that if 7 > 2, then G is growing 
and A > 0. 

The rest of the quantities behave as 

1 a 

H=(ai+a 2 + l)- = -, 
d ( 1\ 1-a 

?= dtUj- i= ^ <o ' Va < 

o 2 = - (af + a 2 + 1 - a\a 2 — a\ — a 2 ) . 



7 



For example we may see in the following table the be- 
haviour of the main quantities for a particular value of 
a 2 . 



a 2 




LO 


7 


G 


A 


1/2 


0.833 


1 


> 2 


/ 


> 


1/2 


0.833 


1/3 


> 2 


/ 


> 


1/2 


0.833 





> 2 


/ 


> 


1/2 


0.833 


-1/3 


< 2 


\ 


< 



choosing other parameters we find another behaviour. 

For this solution we have obtained the following be- 
haviour for the main curvature quantities. We start 
by calculating the curvature invariants, (Jj)^_ , they 
yield: I = 2AT ^ 2 , h = AK x r A , and J 2 = 2K 2 t~ 4 , 
where AT, = K (ax, a 2 , m, n) are numerical constants. 
As it is observed we have obtained a singular solu- 
tion, as it is expect for a SS solution. The non- 
zero components of the Weyl tensor are: Ctxtx = 

Y~ +2(ai — - 1) „— 2mz r~l y +2a\~X„ — 2mz 

Gtyyz — JC/xt 2ai l e 2nZ , Ctztz — 

. . = /C 7 i 2ai e- 2mz , 
gv—e — , where JCi — K,(ax,a 2 ,m,n) 
are numerical constants. Therefore the Weyl invariant 
yields, I3 — (4/3) K^,t~ A . The electric invariant yields: 
I4 = (1/6) AT^" 4 , while the magnetic invariant, 7s, be- 
haves as, I5 = (1/6) K 5 t~ 4 . The Weyl parameter yields, 
W 2 = const., (note that W 2 — > is really small, so the 
model isotropize) while the gravitational entropy yields, 
P 2 = const., as it is expected. 



B. Bianchi III solution. 



This solution is a Bianchi III metric with: 



JC 3 t 2 < ai - 1 >e 2n », 

■■xy 

and C, 



C — r ( . p 2z(n-m)^2(ai+a 2 + l) f< 

^xyxy '^o v , 

— K a t 2a2 e 2nz 
^yzyz — 



ax = ax, a,2 = l, m = 0, n=yl — af, 

so ax £ (0, 1) . In this way the metric collapses to the 
following one: 

ds 2 = -dt 2 + t 2ai dx 2 + t 2 (e 2nz dy 2 + dz 2 ) , (63) 

note that if ax — * 1, then we get 

ds 2 = -dt 2 + t 2 (dx 2 + dy 2 + dz 2 ) , 

this is fact the self-similar solution obtained from a 
Bianchi V metric, which is only valid for lo = —1/3, (in 
fact, the SS Bianchi V, is the FRW solution with k = —1, 
and lo = —1/3.) 

With regard to the main quantities we find that. 

tf = (a 1 + 2)i = |, 

_ d f 1 \ _ 1 - a 

a 2 = 1 -(a 1 -l) 2 . 



Hence we have the following behaviour for the main quan- 
tities: 7 = (w + 1) a £ (1, 6) , a = (ax + 2) £ (2, 3) , and 
A = a\a 2 + ax + a 2 + K > 0, therefore 



p = po^ 



l£ (1,6), 



is a positive decreasing time function. 7 < 2 -4=> to < 0. 
G behaves as: 

G = G i 7 " 2 , 

so it will be a growing function if 7 > 2, constant if 7 = 2 
and a decreasing time function if 7 < 2. 
The cosmological constant behaves as 



A = A 1 



2 

7 

hence its sign will depend on the value of 7. If 7 < 2 then 
we get a negative A, it will behave as a true constant if 
7 = 2 and we find a positive time decreasing function if 
7 > 2. 

For this solution we have obtained the follow- 
ing behaviour for the main curvature quantities. 
We start by calculating the curvature invariants, 
(7i) 2 =0 , they yield: I = 2 (a 2 + ax + 1 - n 2 ) t~ 2 , 
Ix = 4 (3af - 2a? + 2a{ + n 4 - 2n 2 + l) t" 4 , and I 2 = 
2 (a\ + 2a\ + 2a x + n 4 - 2n 2 ai - 2n 2 + l) i~ 4 . As it is 
observed we have obtained a singular solution, as 
it is expect for a SS solution. The non-zero com- 
ponents of the Weyl tensor are the following ones: 



Cf, 



= (1/3)0 



2(01 



a 



tyty 



-(l/6)/Ce 2 



= -(l/6)/C, Cxyxy = (l/Q)JCe 2n H 2a \ Cxzxz = 

(l/6)/Ct 2ai , and C yzyz = -(l/3)0 2 e 2n *, where 
K = (n — 1 + ax) (n + 1 — ax) ■ Therefore the Weyl in- 
variant yields, I3 = (4/3) K, 2 t~ 4 . The electric invariant 
yields: I4 = (l/6)/C 2 i~ 4 , while the magnetic invari- 
ant, Is, vanishes. The Weyl parameter yields, W 2 = 
iC 2 /36(2 + ax) 4 = const., (note that W 2 -> is really 
small, so the model isotropize) while the gravitational 
entropy yields, P 2 = const., as it is expected. 



VII. STRING COSMOLOGICAL MODEL. 

The exponential expansion of the Universe (inflation- 
ary era) causes the Universe to heat up to a very high 
temperature so the subsequent evolution of the Universe 
is exactly as in hot BB model. Hence, the phase transi- 
tion (as the temperature falls below some critical temper- 
ature) in the early universe causes topologically stable de- 
fects: vacuum domain walls, strings and monopoles (see 
[H| and [HI). But domain walls and monopoles are dis- 
astrous for the cosmological models. Strings, on the other 
hand, causes no harm, but can lead to very interesting 
astrophysical consequences (see [35[). Also the existence 
of a large scale network of strings the early universe does 
not contradict the present-day observations. The vac- 
uum strings may generate density fluctuations sufficient 
to explain the galaxy formation (see (36|). 
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The relativistic treatment of strings was initiated by 
Letelier (see [37]]- HH) and Stachel (see [3!|). Here we 
have considered gravitational effects, arisen from strings 
by coupling of stress energy of strings to the gravitational 
field. Letelier (see [38J) defined the massive strings as the 
geometric strings (massless) with particles attached along 
its expansions. Recently, Tikekar and Patel [4l| have 
studied a magnetic string Bianchi type VI model, Bali et 
al [42j have studied a magnetized bulk viscous massive 
string model and Pradhan and Bali [43| have worked with 
a magnetized bulk viscous string model with a variable 
A. 

The energy-momentum tensor, T^ , for a cloud of mas- 
sive strings is given by 



TP 



(p + p) v?m + pgl — Xx^Xi 



(64) 



where p(t) is the rest energy density, i.e. is the rest en- 
ergy density of the cloud of strings with particles attached 
to them (p— strings). A(t) is the string tension density, 
which may be positive or negative, u l is the four-velocity 
for the cloud particles. x l is the four-vector which rep- 
resents the strings direction which is the direction of 
anisotropy and p — p p + A, where p p denotes the par- 
ticle energy density (is the cloud rest energy density), 
i.e. the string tension density is connected to the rest 
energy p for a cloud of strings (p— strings) with particle 
attached to them by this relation (see [44[ and |4a|). 

Since there is no direct evidence of strings in the 
present-day universe, we are in general, interested in con- 
structing models of a Universe that evolves purely from 
the era dominated by either geometric string or massive 
strings and ends up in a particle dominated era with or 
without remnants of strings. 

Moreover the direction of strings satisfy the standard 
relations: 



u 3 Ui = —xPxi = — 1, 
u l Xi = 0, 
x l = (0,a"\0,0) . 



In fact the choice of x l = (0, 0, 0, d _1 (i)) is the only 
one possible since if we choice other direction then from 
the FE we get A = 0. It is customary to assume a relation 
between p and A in accordance with the state equation 
for strings. The simplest one is a proportionality relation 
p = aX, where the most usual choices of the constant 
a are the following ones (see [ilj]). Geometric strings. 
Nambu. a = 1, p = A, p p = 0. p~ strings or Takabayasi 
strings 



a = l + W, 



(66) 



with W > 0, iff p p = WX. Reddy strings. a = —I. A more 
general "barotropic" equation of state is p = p (A) , p p — 
p — X. In the case of Takabayasi strings if W is very small, 
then only geometric strings appear. On the other hand, if 
W is infinitely large, then particles dominate over strings. 



We shall consider the Takabayasi's equation of state, 
i.e. p = aX, with a = 1 + W, so the resulting FE are: 

a'b' a 1 d! d'V K 



a b 
o/_ _ 
a 



d b d 2 
d! b' 



XJI AH 



d! b' n 2 
~d~b ~ d 2 

ad d 2 
a' b' mn 
~a~~b + ~cP 



X (m — n) = 



8ttGp, 


(67) 





(68) 


— 8irGujp, 


(69) 


SnGcop, 


(70) 


-8ttG (up - A) , 


(71) 






4 


(72) 


0. 


(73) 



We have obtained eqs. (|72M73|) from the condition 
divT = 0. As is it observed eq. (|73| bring us to sim- 
plify the FE, since m = n, and then from eq. (|68p we 
get a new constrain, a = b, so the FE collapse to the 
following ones: 



a' d! m 2 „ 
ad a 



a 



a 



m ~ 
I 2 

o 

m 
l 2 
v 



= SnGujp, 



(74) 
(75) 



p > + p(l + u)(2- + -)=X- 



8irG (cop - A) , (76) 
(77) 



So, following the above procedure, and taking into ac- 
count the SS solution, given by eq. (fTT)) . we may find the 
following solution. From ([77]) we get: 



(65) p = Po r p , 0= (l + w)( ai +a 2 + l)- 



1 



= 2, 



1 + W 

and the algebraic system to solve is the following one: 

a\ (ai — 1) + a\ — m 2 = —Au>, (78) 

' ' (79) 
(80) 



2ai (ai - 1) + af + mT = -A w - 



1 + W 



(1+LU) (2dl + l) 

with 



1 



1 + W 



= 2, 



Po = 



A 
8^G' 



A = af + 2di - ra 2 . 



Therefore we may find the following solution: 

I a 2 (w + 1) + 2aru~ 

ai = d2, m=\ = n, (81) 

y w + 1 

w= _ u(2a 1 + l)+2(a 1 -l) ^ ^ 



LU (2dl + 1) + 2d! - 1 
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Discussion: Since W > 0, and a\ > 0, (a\ {u> + 1) + 2a\to > 0) we get: 
I 



to 


m (u>) 


W{uj) 


1 


^af + ai > 0,Voi > 


-^Xit&r-i * °> Va i e (°' V4) =► w e (0, oo) 


1/3 


^/a? + |a x > 0,Vai > 


-il^f >°> Va i e (l/4,5/8)^^6(0,oo) 





ai 


-^^^O.Vaxe (1/2,1] =4ffe (0, oo) 


-1/3 


^/ai (ai - 1) > 0,Vai > 1 


-&E£ ^ °' V °i e f 1 ' 7 / 4 l W e (°> °°) 


-2/3 


•,/ai (ax -4) > 0,Vax > 4 


-§^>0,V ai e (2.7,4] =>\] 



Therefore, by depending on the equation of state for 
the perfect fluid we have a very different behaviour. The 
only nonphysical solution is lu = —2/3. The solution leads 
us to get the following metric: 



ds 2 



with 



-dt 2 



2ai —2mz 



t 2 dz 2 , (83) 



(1 + u) (1 + 2ai) - 



1 + W 



= 13 = 2, 



P = Pot' 



(84) 



In this case we have not obtained the corresponding 
Bianchi III solution since we get only a solution for n — 
rn 7^ 0. As it is observed the metric (|83| has the same 
structure as the metric (|3"5*j) . so it will have the same 
curvature behavior. 



Note that eq. (j73|) is quite strong since it says us that 
if we are working with a Bianchi III metric, where n = 0, 
then the only one (the unique) possibility is m = 0, in 
such a way that the metric collapses to the Bianchi type 
I. This means that there is not massive string Bianchi 
III solution. 



VIII. STRING COSMOLOGICAL MODEL WITH 
TIME VARYING "CONSTANTS". 



Following the above model we shall try to generalize 
it by considering the possible time variation of the "con- 



stants" G and A. In this case the FE are: 

87rGp + A, 
'b' d' 



b' 

b + 


a' dl 
a d 


d'b' 
+ d b 


K 
+ d 2 




fa' 

TO 


5) 


b" 


d" 


d'b' 


2 


~b ' 


f T 


+ ~d~b 


~ I 2 


d" 


a" 


a' d' 


o 

to 


~d 4 


a 


a d 


' ~d? 


b" 


a" 


a' b' 


mn 




a 


a b 


' 1? 




p' 


■HP(1 + 
A (m 


uo)H 
— n) 



A 



(85) 

b dj> ^ 

~8ttGlu P + A, (87) 

SnGujp + A, (88) 

-8vrG {up - A) + A, (89) 

o, 

SnG'p. 



(90) 

(91) 
(92) 



where we have considered the condition div T = 0, eqs. 
(I90H91[) . and therefore we have obtained eq. (|92p as ad- 
ditional restriction. 

As above, we may simplify the FE, so they yield now 
as follows 



a 
a 

a^ + d^ 
a d 

a 





2 

TO 


2-- 

a d 


IP 


d' a' 


2 

TO 


d a 


~ ~d? 




2 




777 


\a) 






d'\ 








^) 




A' 



8nGp + A, 
—SitGup + A, 



(93) 
(94) 



—8itG (uop — A) + A, (95) 
df_ 

7' 

8irG'p. 



X 



(96) 
(97) 



i.e. to = ti, and a = b. By considering the SS solution 
given by eq. (fT7j) we find from eqs. (fM]) and (|6"6"1) 



p = pot"*, 7 = (2oi + 1) (w + 1) 



1 



1 + W 

where we shall assume that po > 0. From eq. 
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obtain: 



A 



9 8irG 
At~ 2 ^Pot- 1 



(99) 



where A — a\ + 2ai — m 2 . 

Now, taking into account eq. ([57)) and eq. (|M|) . algebra 
leads us to obtain 



G = G t 7 ~ 2 , 



Gi) 



A 



(100) 



70 47rpo7' 

and therefore the cosmological "constant" behaves as: 

A = A (l - -J r 2 = A t- 2 . (101) 

With all these results, we find that the algebraic system 
to solve is the following one: 



i\ ~~ 2 



3af - 2a x + mf = A 1 



2 (fj + 1) 

7 

2(td + l) 

7 



(102) 



(1 + W) 1 J ' 

(103) 



and where its solution is: a\ = m = 0, so we arrive to the 
conclusion that there is not SS solution for this model. 



IX. CONCLUSIONS. 

In this paper we have studied some Bianchi types 
VIo&III models under the self-similarity hypothesis. We 
have started by reviewing the "classical" perfect fluid so- 
lution already studied by Collins, Wainwright and Hsu. 
Nevertheless, we have obtained two very restrictive solu- 
tions, which to the best of our knowledge are new. The 
first of them corresponds to a Bianchi type III and is only 
valid for ui G (—1/3, 0) in order to make 02 > 0. The sec- 
ond of them corresponds to Bianchi type VI, and is only 
valid for an equation of state lj € (—1/3, 1). Of course 
we have also obtained the "classical" solution, already 
studied by Collins and several authors. For all these so- 
lutions we have studied their curvature behaviour, i.e. 



we have studied their curvature invariants together with 
the Weyl tensor, the Weyl parameter, W, (showing that 
the models isotropize since (W — > 0)) and we have also 
calculated the gravitational entropy, P 2 . In all the cases 
studied, P 2 = const., since we are working with a SS 
solution. 

In the second of our studied models we have found 
two solutions for the vacuum model. The first of the 
obtained solutions corresponds to Bianchi type VI, and 
to the best of our knowledge it is new. By calculating 
its curvature invariants, (ii)j_ , we have shown that all 
of them vanish, nevertheless the Weyl tensor does not. 
The Weyl parameter is really small and constant and the 
gravitational entropy, P 2 — 00, since I2 = 0. This fact 
tells us that the definition of the gravitational entropy 
does not work for vacuum models. The second of the 
solutions obtained corresponds to Bianchi type III and it 
has already been obtained by Hsu and Wainwright. It is 
a flat solution, for this reason all its curvature invariants 
vanish. 

With regard to a perfect fluid model with time- varying 
constants we have obtained two solutions. Both solutions 
are valid for equation of state u> £ (—1,1]. In this case 
we have been able to enlarge the range of validity for the 
equation of state. The first of the studied solution is of 
Bianchi type VI, while the second one belongs to Bianchi 
type III. In both solutions we have shown that if G be- 
haves as a growing time function then A is a "positive" 
decreasing time function. In the same way, if G is de- 
creasing then A behaves as a "negative" decreasing time 
function. 

For the massive string cosmological model we have 
been able to obtain a solution for Bianchi type VI, where 
the equations of state are ui £ (—2/3, 1] and W £ (0, 00) . 
As we have seen, it is impossible to obtain a Bianchi 
type III solution. Such restriction is obtained from the 
conservation equation, divT = 0. In the last of our stud- 
ied models, which is a massive string cosmological model 
with time-varying constants, we have arrived at the con- 
clusion that there is no SS solution for this model, which 
is quite surprising. 
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